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Abstract 

The inability to predict neutrino masses and the existence of the dark matter are two essential 
shortcomings of the Standard Model. The Higgs Triplet Model provides an elegant resolution of 
neutrino masses via the seesaw mechanism. We show here that introducing vectorlike leptons in 
the model also provides a resolution to the problem of dark matter. We investigate constraints, 
including the invisible decay width of the Higgs boson and the electroweak precision variables, and 
impose restrictions on model parameters. We analyze the effect of the relic density constraint on 
the mass and Yukawa coupling of dark matter. We also calculate the cross sections for indirect and 
direct dark matter detection and show our model predictions for the neutrino and muon fluxes from 
the Sun, and the restrictions they impose on the parameter space. With the addition of vectorlike 
leptons, the model is completely consistent with dark matter constraints, in addition to improving 
electroweak precision and doubly charged mass restrictions, which are rendered consistent with 
present experimental data. 
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I. INTRODUCTION 


The LHC discovery of the Higgs Boson [T] with properties consistent with that of the 
Standard Model (SM) Higgs, while providing a spectacular experimental confirmation of the 
SM, continues to raise questions about SM completeness and about scenarios responsible for 
new physics beyond the SM. In addition, non-collider experimental results confront the SM 
with two major puzzles: neutrino masses and the existence of dark matter. 

The phenomenon of neutrino oscillations shows that at least two neutrinos have nonzero 
but small masses, located around sub-eV scale [2j. The fact that the neutrino flavor structure 
is so different from that of quarks and leptons is a puzzle and may indicate that neutrinos are 
Majorana particles. Many models have been proposed to explain tiny neutrino masses. The 
see-saw mechanism, in which right-handed neutrinos are introduced with large Majorana 
masses M, is perhaps the simplest way to explain tiny neutrino masses. The most direct 
way for implementation of this mechanism for generating neutrino masses is to enlarge 
the particle content of the SM by a complex triplet scalar field, yielding the so-called Higgs 
Triplet Model (HTM) |5j[6]. The neutrino mass problem is resolved at the cost of introducing 
only this additional Higgs representation, together with its associated vacuum expectation 
value (VEV), but without extending the symmetry of the model. The triplet scalar field 
also plays a role in in leptogenesis [7]. 

At the same time, evidence from astrophysics and cosmology indicate that the ordinary 
baryonic matter is not dominant in the Universe. Rather, about 25% of energy density of 
the Universe is comprised of a non-luminous and non-absorbing matter, called dark matter 
(DM). While current observations indicate that most of the matter in the Universe is non- 
baryonic dark matter, they do not provide information on what this dark matter consists 
of. Since the Standard Model (SM), which has been extremely successful in describing all 
current collider data, does not contain any dark matter candidates, a great deal of effort has 
gone into providing viable candidates, or alternatives scenarios (models which include a DM 
candidate naturally). The latter type of models do so at the expense of extra symmetries 
and a much enriched particle content. For models lacking natural candidates, a common 
method is to consider the simplest additions to the SM that can account for the dark matter. 
In these models, the SM particle content is extended by a small number of fields, and a 
new discrete symmetry is introduced to guarantee the stability of the dark matter particle. 
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Several variations can be obtained depending on the number and type of new fields (e.g. 
a scalar, fermion, or vector, a singlet or a doublet under SU(2), etc.) and on the discrete 
symmetry imposed (Z 2 , Z 3 , ...). 

In this work, we look at the Higgs Triplet Model for a resolution to both neutrino masses 
and dark matter problems. The resolution to neutrino masses, alluded to in the above, 
is well-known [5]. The complex triplet couples to left-handed leptons, yielding Majorana 
masses for the neutrinos through L — 2 lepton flavor violating terms |6j, while also con¬ 
tributing to type-II leptogenesis [8j . In addition, extra degrees of freedom that couple to the 
SM Higgs at the tree-level insure cancellation of quadratic divergences to the Higgs mass 
[9], a mechanism where scalars are favored. Additional support for the model comes from 
the observation that heavy particles with strong couplings to the Higgs held can strengthen 
the electroweak phase transition, through the entropy release mechanism from both bosons 
and fermions [10]. 

Unfortunately, as it stands, the Higgs Triplet Model lacks a dark matter candidate. 
Resolutions to this problems were proposed: some with additional Higgs triplets, where the 
neutral component of the additional (real) Higgs representation can act as a DM candidate 
m , another where an additional SU (2 )l triplet scalar fields with hypercharge Y = 1 is 
added [T2] , In this work, we investigate the possibility that the DM candidate is provided 
through the introduction of a complete fourth-generation of vectorlike leptons, comprised 
of SU(2)l doublets plus charged and neutral SU(2)l singlets [13]. A simpler extension of 
the SM with only one fourth generation vectorlikc lepton doublet coupling to a triplet Higgs 
held, which gives Majorana mass to a pseudo-Dirac fourth neutrino has been considered in 

m 

Vectorlike pairs of fermions, unlike their chiral counterparts, are able to have mass ex¬ 
plicitly through the gauge-invariant bilinear interaction in the Lagrangian There 

is no reason why such pairs of vector like fermions do not exist, and many theories such 
as string theories and D-brane theories, often give rise generically to vectorlike states Ha. 
Since the mass of the vectorlike fermions are not generated through the Yukawa couplings, 
the loop contributions involving the Higgs decouple faster than for chiral fermions. Thus 
the constraints from the current Higgs data, precision electroweak observables and direct 
searches are less severe for vectorlike fermions than for chiral fermions. 

Originally, there has been a great deal of interest in vectorlike leptons as a resolution 
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to preliminary data indicating an enhanced Higgs decay rate to diphotons, while the Higgs 
production cross section was in agreement with expectations from SM. The diphoton rate is 
increased through loops of mixed vectorlike leptons. A vectorlike doublet and a vectorlike 
singlet allow for both Yukawa couplings and Dirac masses. The resulting mixing leads to 
a sign flip of the coupling of the lightest lepton to the Higgs field, yielding constructive 
interference with the standard model (SM) amplitude for h —» 77 . 

This does not have to be so in the Higgs Triplet Model, where contributions from vector¬ 
like leptons can be offset by contributions from charged and doubly charged Higgs bosons. 
However new effects of vectorlike leptons can arise. Previous analyses have shown that their 
presence affects the mass bounds and decay patterns of the doubly charged Higgs boson 
116]. improving consistency with present experimental data. 

We extend our previous considerations in [16] to explore the possibility that, introducing a 
new parity symmetry making all new vectorlikc leptons odd, and prohibiting the mixing with 
the ordinary SM leptons, the lightest particle which is odd under this symmetry (a singlet 
neutrino) becomes stable on cosmological timescales, and could have properties consistent 
with it being a candidate for the dark matter of the universe. Note that in a simple heavy 
fourth generation extension of SM, the heavy neutrino does not qualify as a dark matter due 
to its rapid annihilation to SM particles via Z boson exchange na. Leptonic dark matter 
candidates with unsuppressed couplings to the Z boson, such as ordinary fourth generation 
neutrinos are also excluded by limits from direct detection ra- This constraint can be 
relaxed in the model considered here, as the two singlet neutrinos in the model have no 
couplings, or very small couplings, to the Z boson. 

Suppression of the lightest neutrino couplings to the Z boson can also evade present 
experimental limits from LEP on masses of new charged and neutral particles (2|. Measure¬ 
ments of the Z boson width restrict the number of active neutrinos to three, which further 
restrict the mass of the new neutrino to Mn > 39 GeV for a Majorana, and Mn > 45 GeV 
for a Dirac neutrino, precluding the viability of a neutrino which couples to the Z boson as 
a candidate for light dark matter. While, as we will show, we can relax these constraints 
here, the new states will have an effect on the precision electroweak parameters, which we 
calculate and use to restrict the parameter space. We then analyze the consequences of the 
model by requiring consistency with the invisible Higgs width and non-collider experimen¬ 
tal data, particularly with direct and indirect dark matter searches. The relic density, an 
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indication of the abundance of dark matter in the early universe, as measured by PLANCK 
satellite Pj5], is one of the most stringent constraints on any model of DM, as well as direct 
detection experiments search for spin-independent (SI) or spin-dependent (SD) interactions 
with target nuclei, which can be detected by nuclear recoil experiments. Indirect detection 
experiments searches looking for gamma ray excesses measure the annihilation products 
of DM, and their predictions must also be tested in a model of DM. Finally, ultra-high en¬ 
ergy neutrino experiments measure the neutrino flux and flavor composition at astrophysical 
sources. We analyze the predictions for all of these in our model and indicate the constraints 
on vectorlike neutrino mass and coupling which restrict our parameter space. 


Our work is organized as follows. In the next section, Sec. [TT] we summarize the basics 
features of the Higgs Triplet Model with vectorlike leptons. We proceed by examining the 


electroweak precision constraints in the HTM in Sec. Ill, where we present a numerical 
analysis on restrictions coming from the oblique parameters on the masses of the doubly 
charged Higgs bosons and relevant Yukawa coupling. We discuss the invisible decay width 
of the Higgs boson in Sec. IV Then in Sec. [Vj we calculate the dark matter relic density 
and indicate the restrictions it imposes on the mass of the dark matter and on the Yukawa 
couplings. These restrictions are then applied to the evaluation of the spin-dependent and 


spin-independent in the direct detection of dark matter in Sec. VI, and of the annihilation 
cross section of dark matter in Sec. m We discuss detection of DM at colliders in Sec. 


VIII, and then investigate the fluxes of muons and neutrinos from the Sun in Sec. IX We 


summarize our findings and conclude in Sec. |X} 


II. THE HIGGS TRIPLET MODEL WITH VECTORLIKE LEPTONS 


Here we review briefly the HTM with vectorlike leptons, a more detailed version which 
has appeared in our previous work [16], The symmetry group of the HTM is the same as 
that of the SM, with the particle content enriched by (a) the addition of one triplet scalar 
field A with hypercharge Y = 1, and with vacuum expectation value (VEV) va- 


A 


A+ 

V2 

75 (S + VA + irj ) 


A++ 

A+ 

V2 


( 2 . 1 ) 
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and (b) a vectorlikc fourth generation of lepton^J to include: one SU(2)l left-handed lepton 
doublet L' l = {y' L ,e' L ), right-handed charged and neutral lepton singlets, v' R and e' R , and 
the mirror right-handed lepton doublet, L R = (is R , e" R ) and left-handed charged and neutral 
lepton singlets u E and e E , as listed in Table [TJ Note that va is restricted to be small by 
the see-saw mechanism, which requires generation of small neutrino masses, and by the p 
parameter. In general we can assume, conservatively, va H 5 GeV m- 


TABLE I. Representations of vectorlike leptons, together with their quantum numbers under 
SU{3) C X SU(2) l X U (l)y. 


Name 

A 

C"r 

e' r e" l v' r u" l 

Quantum Number 

( 1 , 2 , -1/2) 

( 1 , 2 , -1/2) 

(i, i,-i) (i,i,-i) ( 1 , 1 , 0 ) ( 1 , 1 , 0 ) 


The Lagrangian density for this model contains, in addition to the SM terms, kinetic, 
Yukawa for ordinary leptons, explicit terms for the vectorlike leptons, and potential terms: 

Arm = Ain + A + Al — V (<h, A), (2.2) 


where 


Cy = - + h.c.] 


A AAA A + h.c. 


(2.3) 


are the Yukawa interaction terms for the ordinary leptons, with hfj a 3x3 complex matrix, 
and hij a 3 x 3 complex symmetric Yukawa matrix. Additionally, with the vectorlike family 
of leptons as defined above, 


Al — — \M l L' l L" r + M E e R e" L + M u u R u E + -M' v v' R v' R + -M"v" R v E + h! E (L' L <b)e' R 


+h" E (L" R <S>)e" L + K(L' L r&y R + K' v (L" r t&)v" l + K^L'1%t 2 AL' l + h" L" c zt 2 AL" 
+X E (l' L <S>y R + Ai(Li$)e' i? + Ah LtiT 2 AL' L + A" L^it 2 AL" r + h.c.], (2.4) 


is the Yukawa interaction term for vectorlike leptons and their interactions with ordinary 
leptons, and 


Y($, A) = m 2 $ t $ + M 2 Tr(A f A) + [//$ T iT 2 A f $ + h.c.] + Ai($ t $) 2 

+ A 2 [Tr(A f A)] 2 + A 3 Tr[(A t A) 2 ] + A 4 ($ t <I>)Tr(A t A) + A 5 $ t AA t $, (2.5) 


1 We assume vectorlike quarks to be heavy [Jj and decouple from the spectrum. 
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is the scalar potential for the SM doublet $ (($) = ^|) and triplet A Higgs fields. The 
triplet and doublet Higgs VEVs are related through v 2 = v% + 2v\ ~ (246 GeV) 2 . The scalar 


potential in Eq. (2.5) induces mixing among the physical states for the singly charged, the 


<P 


( 2 . 6 ) 


CP-odd, and the CP-even neutral scalar sectors, which are always small 0(v a/v<&) for the 
first two sectors, but not necessarily so for the latter one, 

( cos a — sin a \ I h 
sin a cos a / \ H 
where the mixing angle is given in terms of the parameters in V (<f>, A) as 
+ o 2n|(A 4 + A 5 ) -4(vIh/V2v a ) 2 

tan 2a = ----- -= ---. (2.7) 

v * 2n|Ai - ( vIh/V2va ) 2 - 2v 2 A (X 2 + A 3 ) 

In our previous work [21], we showed that the Higgs masses and coupling strengths are 
consistent with choosing h to be the SM-like state at 125 GeV, while the state H a lighter 
state, perhaps the state observed at LEP [220 The masses of the neutral h and H are given 
by: 

(rllV V^va) 2 + 2.v 2 A (X 2 + A 3 ) 


mi = 2vi\i cos a + 


m 2 H = 2v\Xi sin 2 a + 


i v $p/V2v A y + 2 u^(A 2 + A 3 


sm a + 


cos 2 a 


va 

- n$n A (A 4 

+ A 5 ) 

sin 2a, 



( 2 . 8 ) 

1 

to 

- v^v A {X 4 

+ A 5 ) 

sin 2a. 

v A 






(2.9) 


The expressions relating the A 1 -A 5 parameters to the Higgs masses can be found in [21J. In 
particular, the doubly charged Higgs boson mass is: 

As q / /i Ar 


2 _ V%H 2 A5 2 _ 

A 3 2 V *~\yft VA 


m H++ = 


5 \ 2 
~2 ^ V,^,, 


( 2 . 10 ) 


where we used v A <C As we choose v A — 1 GeV for consistency with the value of the p 
parameter, the doubly charged mass is approximately m#++ ~ [p — As/v^l (207 GeV). 


The coupling A 5 is expected to be < 1 and for light doubly charged masses, the p parameter 
is smal0 p ~ va ~ Cl (GeV). 

New symmetries can be introduced to restrict the interactions of the vector leptons. For 
instance, we can impose (i) a symmetry under which all the new SU{2) singlet fields are odd, 


This scenario was imposed by the requirement of an enhanced diphoton signal for the Higgs of mass 125 


GeV, so it can be relaxed here. 

3 Specifically, for our parameter space p = 0.2 GeV and A 5 < 0. 
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while the new SU(2) doublets are even, which forces all Yukawa couplings involving new 
leptons to vanish, h' E = h" E = h' v = h" = hX = hX = 0, and the vector lepton masses arise 
only from explicit terms in the Lagrangian [13] : and/or (ii) impose a new parity symmetry 
which disallows mixing between the ordinary leptons and the new lepton fields, under which 
all the mirror fields are odd, while the others are even |23], such that X E = X l L = XX = XX = 
A" = 0. The latter are important for light vectorlike leptons, as this scenario would satisfy 
restrictions from lepton-flavor violating decays, which otherwise would either force the new 
leptons to be very heavy, ~ 10 — 100 TeV, or reduce the branching ratio for the Higgs into 
dileptons to 30-40% of the SM prediction. In addition, if all vectorlike leptons are odd under 
this symmetry, the lightest particle can become stable and act as all, or part of, the dark 
matter in the universe. Thus the assumption (ii) has all the attractive features we like for 
this analysis, and we adopt it here, while allowing h' E , h E , h ' y , h", hX, hX % 0. 

As we concentrate on the possibility that the lightest neutral component of the new 
vectorlike leptons is a dark matter candidate, we are primarily interested in light states. 
The 2 x 2 mass matrix M e for the charged sector is defined as using 


(e' L %)(M e ) 



with 


■Me — 


I m' E M l 
\M e m" E 


( 2 . 11 ) 


with m' E = h' E v®/y/2 and mI E = h" E v<&/y/\ 2, from the Lagrangian Eq. 
can be diagonalized by two unitary matrices U L and U R as follows: 


(2.4). 


The mass matrix 


u l] m e u r 


( Me ! 0 

y 0 Me 2 


( 2 . 12 ) 


The mass eigenvalues are (by convention the order is M El > M E2 ) 


m e u e 2 


(M, 


l + < + Ml + m" E 2 


)± 


/ 2 


(Mi + m‘ E 


- Ml - m" E ) 


+ 4 {m" E M L + m' E M E ) 
(2.13) 


while in the neutral sector the mass matrix is: 


1 

2 



U R V R C 


( U ’ lC \ 

v'r 

\< C ) 


with M v 


1 0 

m' v 

M l 

0 ^ 

™>'v 

M' 

■LV-Ljy 

0 

My 

M l 

0 

0 

m" 

\ o 

M v 

m" 

k) 


<d) (Vl„) 


(2.14) 










with m' v = h' v v*/y/2 and m" = h"v$/y/ 2. This mass matrix can be diagonalized by a unitary 
matrix V: 


V f M„V 


(M U1 0 0 

0 M V2 0 


0 ^ 
0 


(2.15) 


0 0 M U3 0 

o 0 0 M V J 

In the limit where the explicit mass terms Ml, M R and M v in the interaction Lagrangian 
vanish, after electroweak symmetry breaking there are two charged leptons with masses m' E 
and 77i e, and four neutrinos with masses: 





+ m ! 2 ± 




+ m’l 2 ± 



(2.16) 

(2.17) 


The lightest of these eigenvalues will be the dark matter candidate, and as it is odd under 
the additional parity symmetry (n), it is stable. For vanishing h' u , h'l Yukawa couplings, 
the two singlet vectorlike neutrinos have vanishing couplings with the Z boson. Lifting 
the Yukawa couplings slightly from 0 allows mixing between the singlet neutrinos and the 
neutral components of the doublet vectorlike leptons, inducing a (small) coupling to the Z 
boson. For simplicity, we adopt the scenario in [T3] where h' v ^ 0, but h'l = 0, as well as 
setting the explicit neutrino mass in the Lagrangian M u = 0. This scenario is sufficient 
to provide a single DM candidate and a single Yukawa coupling, and transparent enough 
to yield consequences. It corresponds to one neutrino state which does not mix and is 
sterile (the mirror SU(2)l doublet u R ), while the remaining neutral sector consists of three 
neutrinos which mix, with mixing matrix in the {v'l, u' R , u R ) basis given by 


^ 0 m' M l ^ 




ml Ml 0 
M l 0 0 


(2.18) 


The Yukawa coupling h' u must remain small to insure smallness of couplings to the Z boson. 
In the limit h' v = 0, the matrix has two degenerate eigenvalues of mass Ml, predominantly 
SU(2)l doublets, and one state with mass M' and predominantly singlet. For h! v ^ 0, 
these three states mix, generating a small mixing coupling to the Z boson. The lightest 
neutrino state M U1 emerges as being dominantly v' R and is the dark matter candidate. For 
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the charged lepton sector, we take Ml = 205 GeV and Me = 300 GeV and h' E = h" E = 0.8 
[Qj. In this case the lightest charged lepton will be M E2 ~ 108 GeV, close to the LEP 
limit, M e > 102.6 GeV [2], which imposes an upper limit on the mass of the dark matter 
candidate Mdm = M Vl < Me 2 - 

Next, we analyze the effects of the new states on electroweak precision parameters in the 
HTM and consequently, the restrictions imposed on its parameter space. 

III. VECTORLIKE LEPTON CONTRIBUTIONS TO THE S AND T PARAME¬ 
TERS 

Adding new particles to the model spectrum affects quantum corrections on the propa¬ 
gators of W and Z bosons. The corrections are parametrized by two oblique parameters, S 
and T [j] which encapsulate the model restrictions coming from electroweak precision data. 
For a Higgs state with mass rrih = 125 GeV, the allowed ranges are [2D] 


AS — S — Ssm — 0.05 ± 0.09, 
AT = T-T sm = 0.08 ± 0.07, 


(3.1) 


with errors correlated by a factor of 0.88. The explicit expressions for the S, T and U 
parameters for the HTM are given in [ID]. The addition of vectorlike leptons modifies these 
by the following contributions. For the S parameter ra. 


7r 


1 



3 


3 



(3.2) 


4 We set U = 0. 
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while the oblique correction parameter T for vectorlikc leptons is [13] : 


2 3 


T = 


47 r S W C W M Z 


- 2 EE< i t/ triu t i 2 +i^rin l i 2 )6 3 (M„ t ,M Ej ,o) 

j=1 k= 1 


2 3 


2 E E0) 

j=l fc=l 

E (inf iwr + infim 2 ) oj 


3 


j,fc=i 

+ (l^nl 4 + + (If/^l 4 + 


+ (2|^iil1^ir + 2|C/"nt/"| 2 )6 3 (M El ,M E2 ,0) 

- £ ReiU^U^U^M^MEMME^M^O) \ , 
j,fc=i 


(3.3) 


where the Passarino-Veltmann functions are: 


b 0 (M 1 ,M 2 ,q 2 ) = 

f 1 A 

/ 

(3.4) 

bi(Mi, M 2 , q 2 ) — 

r 1 a 

J xl og( — )dx, 

(3.5) 

b 2 (M 1 ,M 2 ,q 2 ) = 

J x(l-x) log(^)dx = b 2 (M 2 ,M 1 ,q 2 ), 

(3.6) 

6 3 (Mi, M 2 , 0 ) = 

M 2 , 0) + M 2 b\ (M 2 , M u 0) 

(3.7) 

2 

h (m 1 ,m 2 ) = 

M 4_ M 4 + 2M 2 M 2 log(^i) 

MlM2 2(M 2 - M|)3 

(3.8) 


We defined A = M%x + Mf( 1 — a;) — x(l — x)g 2 and in the above A 2 is an arbitrary 
regularization scale that will not affect physical observables. The function , M 2 ) = 

— 1/6 remains well defined in the limit M 2 —> Mx. As in the HTM without vectorlike 
leptons, the S parameter does not impose any restrictions on the parameter space of the 
model, while the T parameter is very restrictive. The reason is that T depends quadratically 
on mass differences, while S only logarithmically. 

We proceed to analyze restrictions on the relevant masses and couplings in the model 
coming from the T parameter. In Fig. [l] we show the effects on the T parameter as a 
contour in a m#±± — sin a plane (with since the mixing angle in the neutral Higgs sector), 
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Sin a 


for two values of dark matter masses, Mom = 30 GeV and Mom = 50 GeV. The allowed 
values for this parameter, —0.2 < AT < 0.4, are given in the code bars (colored contours in 
the figure). The maximum doubly charged mass values allowed are m H ±± ~ (280 — 290) GeV 
for dark matter masses in 24-30 GeV and 70-90 GeV regions, and approximately 250-270 
GeV for dark matter masses in 30-70 and 90-103 GeV range. We have selected the particular 
values for Mom = 30 GeV and h' u = 0.65 (left panel) and Mom = 50 GeV and h' u = 0.28 
(right panel) to belong to the parameter space where the relic density is within experimental 
bounds, as explained in detail in Sec. |Vj As the figure indicates, the T parameter depends 
only slightly on sin a, but is extremely sensitive to the mass of the doubly charged Higgs 

0.5 

0.0 

-0.5 

- 1.0 

0 50 100 150 200 250 300 0 50 100 150 200 250 300 

M oi ++ M oi ++ 




FIG. 1. (color online). Contour graphs showing the contribution to the T parameter in the HTM 
(as given in the code bars) with vectorlike leptons, as a function of the doubly charged Higgs mass 
mjj±± and the mixing angle sin a, for fixed values of the neutrino Yukawa coupling h' v . We take 
(left panel) Mom = 30 GeV, h' v = 0.65, (right panel) Mom = 50 GeV, h' u = 0.28. The allowed 
range of T parameter is —0.2 < AT < 0.4. The white region represents the parameter region ruled 
out by the constraints. 


In Fig. [2] we show the variation of the T parameter as a contour in an rn H ±± — h' v 
plane (left panel) and h' v — sin a (middle and right panels). For the left panel, we chose an 
illustrative example with Mom = 30 GeV and sin a = 0.5; the T parameter does not depend 
sensitively on varying these, but again it is very sensitive to the mass of the doubly-charged 
Higgs boson, as shown in the middle and right side panels, where increasing the value of 
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FIG. 2. (color online). Contour graphs showing the contribution to the T parameter in the HTM 
with vectorlike leptons (values as given in the code bars, within the allowed range —0.2 < AT < 
0.4.) as a function of the parameters of the model. In the left plot, we show the combined 
dependence on h' v and m H ±± (for Mdm = 30 GeV and sin a = 0.5); in the middle (right) panel, 
the dependence on h' v and sin a for Mjjm = 50 GeV and mjj±± = 240 GeV (mjj±± = 260 GeV). 
The white region represents the parameter region ruled out by the constraints. 


mjj±± from 240 to 260 GeV places significant restrictions on the T parameter. Increasing 
the mass of the doubly charged Higgs boson and decreasing sin a (the mixing angle) impose 
restrictions on h! v (the vectorlike neutrino Yukawa coupling) from the T parameter. Note 
that the T parameter is not sensitive to the mass of the dark matter candidate and it affects 
it only indirectly, through the restrictions on the Yukawa couplings. 


IV. INVISIBLE DECAY WIDTH OF THE HIGGS BOSON 

The existence of the vectorlike neutrino v\ as dark matter candidate will have an effect 
on the branching ratio of the Higgs boson, if > 2 M U1 . Given that Ui is stable, the 
decays h —> h — > UiPi will contribute to the invisible Higgs branching ratio, which 

is constrained by combined CMS and ATLAS measurements to be BR my < 58% for a SM 
Higgs with a mass of 125 GeV m and more stringently by global fits to be BR nrv of 29% 
with 95% C.L. |25j . 

In the Higgs Triplet Model, the tree-level decay width of the Higgs boson into vectorlike 
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neutrinos is [2TJ [26] 


where 


[r(/r —* iq^i)]ffTM — 


G F m h {M Pl C ^ ll71 ) 2 / _ 2 M u 2 \ 2 
2tTa/ 2 V wifc / 


cos 2 a, 


= ^KRe(V u Vn) 


(4.1) 

(4.2) 


is the Higgs coupling to the lightest vectorlike neutrino (iq). As well, the component from 
the neutral triplet Higgs field violates lepton number and can decay into two neutrinos as 


[r (h —>• viVi)}htm = r [h —>■ v[ui) + T(/r —> t'fk'i) 


= V 

2 4vr I ml 


Ml 
1*4 —^ 
ml 


sm a, 


(4.3) 


where /i(, is the triplet coupling constant from Eq. (|2.4|). The invisible branching ratio of 


the Higgs boson is defined as 


BR mv = 


[T(h —> viVi )}htm + [r(/i —> viVi)\htm 


(4.4) 


[T(h —>• ^i^i)]nTM + [r(h —> zqi'OJ.HTM + [r(h)]/fTiu'’ 
where [r(/r)]^ TM is the total Higgs decay width in the HTM without vector like leptons. 

In Fig. [3j we show the invisible branching ratio of Higgs boson ( BR imr ) in the HTM 


with vectorlikc leptons as a contour plot in an M dm — h' u plane, for triplet Yukawa coupling 
h! i = 0.01. We compare the calculation with the upper limit on BR my derived from 
global fits to ATLAS and CMS data |25|^) As expected, the region restricted is only for 


Mdm < Wh/V, where the Higgs can decay to pairs of dark matter with a sizeable width. 
The left panel depicts the invisible width for the mixing angle in the neutral CP-even Higgs 


sector, sin a = 0.1, the middle panel for sin a = 0.5 and the right panel for sin a = 0.8. The 
figures show that increasing the Yukawa coupling (h' u ) results in an increase the invisible 
branching ratio of Higgs boson ( BR irw ) as the decay into DM is enhanced, while decreasing 
sina- imposes more restrictions on h' u in order to get correct BR irw , indicating that both the 
doublet and triplet Higgs components play an important role in the invisible decay. 


V. DARK MATTER RELIC DENSITY 

Global fits to a number of cosmological data (Cosmic Microwave Background, Large Scale 
Structure and Type la Supernovae) determine very precisely the amount of non-baryonic 
5 These global fits, though more restrictive, are completely consistent with our analyses and do not restrict 
the parameter space unnecessarily. 
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FIG. 3. (color online). Contour graphs showing the invisible branching ratio of Higgs boson 
(BRi nv ) in the HTM with vectorlike leptons, as functions of the dark matter mass Mom = M V1 
(GeV) and the neutrino Yukawa coupling h! v , for h! v v = 0.01 . We compare to the upper limit of 
BR[ nv = 29 % from global fits to ATLAS and CMS data [251 and we chose sin a = 0.1 (left panel), 
sin a = 0.5 (middle panel), sin a = 0.8 (right panel). If the mass of DM neutrino is in the white 
region, it does not contribute to the Higgs invisible decay width. 


(DM) in the energy-matter of the universe at floM^ 2 = 0.1123 ± 0.0035 [27] . where Qdm is 
the energy density of the DM with respect to the critical energy density of the universe, and 
h is the reduced Hubble parameter. Any analysis of DM must correctly replicate this value. 

To this end, we used CalcHEP [28] to implement the Lagrangian of the HTM with vec¬ 
torlike leptons into micrOMEGAs [22], to calculate the relic density (CLom^ 2 ), spin-dependent 
cross section ( a SD ), spin-independent cross section ( a SI ), annihilation cross section ((cry)), 
and the flux of neutrino and muon predicted by the model. For the purpose of comparing 
with the data, we consider the 2cr allowed range of relic density: 0.1144 < TIom^ 2 < 0.1252, 
as constrained by WMAP [27] and PLANCK pL9| . 

In Fig. [4] we present the allowed range of relic density of dark matter as function of the 
dark matter mass M dm (GeV) and the Yukawa coupling h' u , for two different values of the 
mixing angle, sin a = 0 (left panel) and since = 0.8 (right panel). Due to resonant annihila¬ 
tion into Z bosons or Higgs boson h respectively, we can see two dips at Mom ~ 45 GeV 
and Mom ~ 62 GeV. For a fixed Yukawa coupling (h() the cross sections becomes enhanced 
at the Z pole and similarly at the Higgs pole, with a dominant decay into quark/antiquark. 
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As the dark matter relic density is inversely proportional to the annihilation cross section, 
the relic density decreases in these regions. Thus, in order to produce the correct dark 
matter relic density, we need to decrease the value of Yukawa coupling h' v to compensate 
for the effects of Z and h resonances, which produces the two dips at Mdm = Mz /2 and 
Mdm = Trih/2. Above Mdm = 80 GeV annihilation into W + W~ pairs (and later also Z 
bosons) becomes kinematically accessible. Finally, the relic density becomes dramatically 
suppressed for M dm ~ 100 GeV due to co-annihilation with the lightest charged vector¬ 
like lepton [El EBj- The effect of the Higgs resonance at M dm ~ 62 GeV is slightly more 
pronounced for sin a = 0.8 than for sin a = 0 (this is the effect of increasing the triplet 
component contribution) and, above Mdm = 80 GeV, for sin a = 0.8, the relic abundance 
decrease is slightly more pronounced than in the case with sin a = 0, but the changes are 
small. Overall, the graph for sin a = 0.8 shows no marked difference from the one with 
sin a = 0. The results shown are for m#±± = 240 GeV. We calculated relic density for 
different values of the doubly-charged Higgs boson mass and found that it is insensitive to 
variations in this parameter. Relic density constraints restrict the dark matter mass to be 
heavier than 23 GeV and lighter than 103 GeV in our model, independent of any other 
parameters, such as Yukawa couplings or mixing angles. 


VI. DIRECT DETECTION 

Dark Matter is spread over the whole universe. This provides the opportunity to detect 
it as it passes through and scatters off normal matter (neutrons or protons), producing 
detectable signals. Though direct detection is the most straightforward method of detecting 
DM, such events are very rare, deposited energies very small, and thus direct detection 
requires very sensitive detectors with highly accurate background rejection. The expected 
signals depend on the nature of the DM. For vectorlike neutrinos, annihilation through 
the Higgs or Z boson exchange is expected to yield significant rates for direct detection. 
The interaction of DM with nuclear matter can be classified as elastic or inelastic; and as 
spin-dependent or spin-independent. 

In elastic scattering the DM interacts with the nucleus as a whole, causing the nucleus 
to recoil, while in inelastic scattering some of the energy goes into recoil, and some is used 
to excite the nucleus to a higher energy state, from where it decays by emitting a photon. 
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FIG. 4. (color online). Contour graphs showing the correct relic density of dark matter as function 
of the dark matter mass Mdm (l n GeV) and the neutrino Yukawa coupling h! v in the HTM with 
vectorlike leptons for sin a = 0 (left panel) and sin a = 0.8 (right panel). We impose the restriction 
0.1144 < Qouh 2 < 0.1252. The relic density is insensitive to the doubly charged Higgs boson 
mass, chosen here to be 240 GeV. 

The dark matter detection experiments (DAMA/LIBRA [30], CoGeNT [21] and CRESST-II 
02!) have reported signals consistent with a light DM candidate and with an elastic cross 
section with nucleons of (9(10 -41 — 10 _4o cm 2 ). 

In spin-dependent (axial vector) scattering, the DM spin couples with the spin of the 
nucleon, while in spin-independent (scalar) scattering, the cross section does not depend on 
this, and thus it is larger for larger nuclei because of the coherence of DM interacting with the 
nucleus as a whole. We analyze the predictions of our model for the spin-dependent and spin- 
independent cross sections in turn, and compare them with the experimental predictions. 

In Fig. [5j in the upper panels, we present the spin-dependent (SD) cross section of dark 
matter scattering off nucleons, as a function of the dark matter mass M dm for sin a = 0. 
The left panel is for the proton, the right one for the neutron. The red lines show points 
of the parameter space, with restricted Mdm and h! v values, which reproduce acceptable 
relic density. The areas above the pink dashed line and green dashed-dotted line are ruled 
out by the COUPP [33] and XENONIOO [34] measurements, respectively. As the plots 
show, to obtain the correct relic density, the resonantly enhanced annihilation rate implies a 
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suppressed Yukawa coupling for the neutrino DM, which leads to a suppressed cross section. 
Here again we observe the two dips surrounding the Z resonance and the h resonance. 
The limits on the SD cross section from COUPP and XENONIOO results do not restrict 
the parameter space of our model. In the bottom panel, we plot contour graphs for the 
spin-dependent cross sections of the nucleon as functions of the dark matter mass Mdm 
and Yukawa coupling h' u , for sin a = 0. Again we show the spin-dependent cross section of 
the proton and neutron in the left and right panel, respectively. All points are consistent 
with experimental bounds on the spin-dependent nucleon cross sections, as indicated by 
the color-coded panels, but only parameter points situated along the dashed-dotted lines in 
the bottom panels give the correct dark matter relic density. These cross sections are not 
sensitive to variations in sin a. 

In Fig. |bj we plot the spin-independent (SI) cross section of nucleon, as a function of the 
dark matter mass Mdm (in GeV) for sin a = 0 (left panel). The red line includes all points 
yielding consistent relic density. The regions above dash-dotted black line, dash-dotted 
green line, dash-dotted orange line, dash-dotted blue line, dash-dotted purple line, dash- 
dotted pink line are ruled out by XENONIOO |33], XENONIOO with 2cr expected sensitivity, 
CRESST-II [32], CDMS-II [SB], TEXONO [37] and DAMIC100 (expected for 2014) [38] 
results, respectively. The cross section is enhanced at the Z pole and h pole and there, for 
a suppressed direct rate, the Yukawa coupling must be suppressed to compensate for the 
resonant production effect. This is seen as two dips at Mdm ~ Mz /2 and Mdm ~ m^/2. 
The limit on the SI cross section from XENONIOO constrains strongly our model, while the 
updated results from the other experimental results do not restrict the parameter space. As 
the left panel of the figure shows, XENONIOO results (with 2cr expected sensitivity) restrict 
the dark matter mass to be in the 37-52 GeV, or 57-63 GeV ranges, or heavier than 95 GeV. 
In the middle panel, we show the spin-independent cross section of the proton as a graph 
in Mdm ~ h' u space, constrained by all the experiments with the exception of XENONIOO, 
while in the right panel we include XENONIOO measurements. The latter rules out large 
regions of parameter space (in white) while in both panels colored contours (as coded in the 
attached bars) are allowed by the spin-dependent experiments. In both middle and right 
panels, the dashed-dotted line represents the only parameter points with acceptable relic 
density. Note here that, in agreement with the left panel, there are regions of the parameter 
space where no combination of Mdm and h' v satisfy both relic density and XENONIOO SI 
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FIG. 5. (color online). Top: the spin-dependent cross section of the nucleon in the HTM with 
vector-like leptons, as a function of the dark matter mass M^m (GeV) for sin a = 0. We show 
(left panel) the spin-dependent cross section of the proton (red line) with XENON100 (dash-dotted 
green) and COUPP (dashed pink) io4\/ results, (right panel) the spin-dependent cross section of 
the neutron (red line) with XENON 100 (dash-dotted green) results (34J. The area above the 
pink dashed line and green dash-dotted line are ruled out by the COUPP and XENONIOO results 
respectively. Bottom: Contour plots showing the spin-dependent cross section of the nucleons in 
the HTM with vectorlike leptons, as functions of the dark matter mass Mjom and Yukawa coupling 
h! v , for sin a = 0. We show the spin dependent cross section of the proton (left panel) and the 
spin-dependent, cross section of the neutron ( right panel). The panels at the right indicate the 
color-coded values of the cross section along each slice, and the dashed-dotted line represents the 
only parameter points with acceptable relic density, 0.1144 < < 0.1252. 
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cross section restrictions. Here too, the cross sections are not sensitive to the mixing angle 
or to other parameters in the model, and for the spin-independent, the cross sections for the 
proton and neutron are indistinguishable. 
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FIG. 6. (color online). Left panel: the spin-independent cross section of the proton as a function 
of the dark matter mass Mum (GeV) in the HTM (red line). We also show XENONIOO j35j / 
(dash-dotted black) XENONIOO with 2o expected sensitivity (dash-dotted green), CRESST-II 
J32j / (dash-dotted orange), CDMS-II 136\j (dash-dotted blue), TEXONO j3ff (dash-dotted purple) 
and DAMIC100 (expected in 2014) 13#] / (dash-dotted pink) results. Middle panel: Contour graph 
showing the spin independent cross sections of nucleon in the HTM with vector-like leptons, as 
functions of the dark matter mass Mum and h' v for for sin a = 0, considering all experimental 
constraints except XENONIOO. Right panel: same as the middle panel, but including constraints 
for XENONIOO (with +2o expected sensitivity) |35f upper limit. The values of the cross section are 
indicated on the plots. The contours indicate points consistent with the respective experimental 
constraints, while the dashed-dotted line includes only points with acceptable relic density, 0.1144 < 
0 DM h 2 < 0.1252. 


VII. INDIRECT DETECTION 

Pairs of dark matter particles annihilate producing high-energy particles (antimatter, 
neutrinos or photons). Indirect detection experiments for dark matter look for signatures 
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of annihilations of DM originating from particles in the flux of cosmic rays and are sen¬ 
sitive to DM interaction with all the SM particles. The most stringent constraints on 
DM annihilation cross sections have been derived from the Fermi Gamma Ray Space tele¬ 
scope (Fermi-LAT) [39], used to search for DM annihilation products from dwarf spheroidal 
galaxies and the Galactic Center, which probe annihilation cross sections into photons of 
(cry) ~ 3 x 10 _26 cm 3 /s. These searches have attracted a lot of attention due to the unex¬ 
pected high flux of cosmic ray positrons observed by the PAMELA experiment [40], and 
confirmed by AMS | [T . 


In Fig. [7J we present the annihilation cross section of DM as a function of the dark matter 
mass Mdm and compare it with the constraints on the dark matter annihilation cross section 
for the e + e~ channel, channel, r + r~ channel, uu , bb channel and W + W~ channel at 

95% CL, derived from a combined analysis of 15 dwarf spheroidal galaxies of Fermi-LAT 
Collaboration results [39] (left panel). As the figure shows, the limit on the annihilation 
cross section from Fermi-LAT Collaboration results does not impose any restriction on our 
model parameters. Again, the annihilation cross section is enhanced at the Z pole around 
Mdm = Mz/ 2. The regions around Mdm = Mz /2 can be brought into agreement with 
the relic density constraint by modifying the neutral Yukawa coupling h' v . In order to treat 
the regions where the annihilation cross section is enhanced, we need to decrease the value 
of Yukawa coupling. A suppressed coupling leads to suppression of the annihilation rates 
Hal EH]. The effect of the Higgs pole at Mdm ~ h/2 is more dramatic than the effect of 
the Z pole. The dominant annihilation modes of dark matter pair in this region are coming 
from decay into quark/antiquark (mainly b b, which gives a relative contribution of ~ 77% 
to 1/VLdm^ 2 ) and also small contribution from cc and rr to obtain the correct dark matter 
relic density. In the right panel, we show the annihilation cross section as a contour plot in 
the dark matter mass Mdm and Yukawa coupling h' u plane. Here, as in the previous figures, 
the contours, according to the color-coding in the attached bar, represent the regions of 
the parameter space consistent with the experimental results, while the white regions are 
excluded. Only points along the dashed-dotted line have acceptable relic density. 
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FIG. 7. (color online). Left panel: the annihilation cross section of DM as a function of the dark 
matter mass Mjjm (GeV) (red line). We show the constraints on the dark matter annihilation 
cross section for the e + e~ channel (dash-dotted green), p + channel (dash-dotted orange), t + t~ 
channel (dash-dotted pink), uu (dash-dotted blue), bb channel (dash-dotted purple) and W + W~ 
channel (dash-dotted black) at 95% CL derived from the combined analysis from Fermi-LAT Col¬ 
laboration m Right panel: contour plot showing the annihilation cross section as function of 
the dark matter mass Mjjm an d its Yukawa coupling h' v . The contours are consisted with the 
experimental values for the cross sections, indicted by the color-coded bar, while the white regions 
are ruled out. Only points along the dashed-dotted line give the correct dark matter relic density, 
0.1144 < D dm h 2 < 0.1252. 

VIII. DETECTION AT PARTICLE COLLIDERS 

If dark matter has significant coupling to nuclear matter it can be produced in high en¬ 
ergy collisions at LHC or at future colliders. Once produced, as it is neutral and weakly 
interacting, DM will not be observed directly, but it could inferred from missing transverse 
momentum. Collider searches provide the opportunity to study DM production in a con¬ 
trolled environment. They are particularly sensitive to the region of low mass dark matter, 
where backgrounds are smaller. At the LHC, dark matter can be produced directly, together 
with additional radiation from the quarks or gluons participating in the reaction, which re¬ 
sults in a single jet (mono-jet) plus missing momentum. High energy lepton colliders could 
create dark matter through a similar process. Assuming DM couples to quarks and gluons 
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and couplings the order of the electroweak size, LHC excludes DM masses up to 500 GeV 
and for DM coupling to electrons with the same-size couplings, LEP excludes DM with mass 
below 90 GeV. Neither of these restrictions are applicable here, as vectorlike neutrinos do 
not couple directly to either quarks or leptons. 

IX. THE FLUX OF MUONS AND NEUTRINOS FROM THE SUN 

The recent observation of ultra-high energy neutrino events at IceCube H23 seem to in¬ 
dicate a possible deficit in the muon track (known as the muon deficit problem) and an 
apparent energy gap in the 3-year high energy neutrino data, challenging a simple expla¬ 
nation in terms of atmospheric neutrinos and suggesting an extra-terrestrial origin. These 
astrophysical neutrinos are assumed to have originated from the decays of charged particles 
produced in pp or py collisions. While the data obtained is largely consistent with SM pre¬ 
dictions, the flux shows a mild deficiency in muons at high energies, prompting alternative 
explanations involving dark matter. 

In Fig. [8j we show the neutrino (left panel) and muon (right panel) fluxes as functions 
of the dark matter mass M dm (GeV). On the top graphs, we plot our results as a red 
curve, and include the upper limits on the neutrino and muon flux for the bb channel, r + T~ 
channel, and the z/ e z7 e , ryzy,, v T v T channels from the Baikal NT200 detector results [43]. 
While the limit on the muon flux from Bailkal results does not impose restrictions on our 
modeQ the neutrino flux excludes DM particles with mass in the 74-85 GeV. The figures 
show again the two dips at Mdm ~ 45 GeV and Mdm ~ 62 GeV. Unlike in the annihilation 
cross section, here the effect of the Z pole is more dramatic than that of the Higgs pole at 
Mdm h/2. The bottom panels show the fluxes of neutrino and muon as contour plots 
in the dark matter mass Mdm and the Yukawa coupling h' v . Note that here, as before the 
contours are consistent with the experimental values for the measured flux of muons and 
neutrinos. However only points along the dashed-dotted line are consistent with the dark 
matter relic density exclusion limit. 


6 Due to limited space on the figure, we show the recent results of Baikal NT200. But our results are also 
consistent with those from the Baksan Neutrino Observatory [44] . 
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FIG. 8. (color online). Top: the flux of neutrinos (left panel) and muons (right panel) in the 
HTM with vectorlike leptons, as functions of the dark matter mass Mjjm (GeV) (red line). We 
also show the upper limits on the neutrino and muon flux for the bb channel (dash-dotted green), 
t + t~ channel (dash-dotted orange), v e v e channel (dash-dotted pink), channel (dash-dotted 
blue) and v T v T channel (dash-dotted purple) from the Baikal NT200 detector results '(T5I . Bottom: 
contour plots showing the flux of neutrinos (left panel) and muons (right panel) in the HTM with 
vectorlike leptons, as functions of the dark matter mass and Yukawa coupling h' u . The flux 

values are color-coded as in the bar attached. The contours indicate points consistent with the 
respective experimental constraints, while the dashed-dotted line includes the only points with 
acceptable relic density, 0.1144 < < 0.1252. 
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X. SUMMARY AND CONCLUSION 


In this work we analyzed the effects of introducing vectorlike leptons in the Higgs Triplet 
Model. Our aim was to provide a scenario that can explain both neutrino masses and provide 
a DM candidate, problems unresolved in the SM. We choose a full generation of vectorlikc 
leptons (one left-handed doublet, two right-handed singlets, together with their mirror rep¬ 
resentations). We insure that a new symmetry differentiates between ordinary leptons and 
the new states, forbidding unwanted lepton flavor violation. Opting for a simplified Yukawa 
coupling structure, a mostly singlet right-handed vectorlikc neutrino emerges as a single 
DM candidate. Introducing vectorlike leptons in the HTM relaxes the severe constraints on 
the mass of the doubly charged Higgs boson coming from electroweak precision tests. We 
revisit precision observables in this work, and show that while the S parameter does not 
impose constraints on the parameter space, the T parameter is restrictive, allowing only 
certain combinations of doubly-charged mass, Yukawa couplings and mixing angles between 
the neutral Higgs bosons. Of these, the most sensitive parameter is the mass of the doubly 
charged Higgs boson, required to be less than about 280 GeV, but this boson has different 
branching ratios than in the minimal HTM. The T parameter is insensitive to the mass of 
the dark matter candidate. 

We verify that the invisible decay width of the Higgs boson is consistent not only with 
the experimental data, but with the more restrictive limits imposed by global fits to the 
Higgs data. The invisible width is a relevant constraint for dark matter masses less than 
1/2 the Higgs mass, and all of these survive. More stringent constraints come from direct 
detection experiments, especially from restriction on spin-independent nucleon cross section, 
and from the relic density. The latter restricts the combination between dark matter mass 
and its Yukawa coupling to narrow bands in the parameter space, and disallows entirely 
regions where the DM candidate is lighter than 23 GeV, or heavier than 103 GeV. If one 
includes constraints from XENONIOO on spin-independent scattering of dark matter off 
nucleons, these further restrict the dark matter mass to be in the the ranges: 37-52 GeV, 
or 57-63 GeV, or heavier than 95 GeV, all for points satisfying relic density constraints. 
In addition, consistent with direct detection experiments, the neutrino flux excludes DM 
particles with mass in the 74-85 GeV range. These are the most stringent restrictions, and 
they are insensitive to other model parameters, such as other masses (particularly the doubly 
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charged Higgs boson) and the mixing angle between the neutral Higgs bosons. 

To summarize, we have presented a simple model that accounts for neutrino masses and 
dark matter and is consistent with the relic density and all direct and indirect searches. This 
model assumes a single dark matter particle, and the experimental data restricts its mass to 
be confined to limited regions in the parameter space. If the dark matter is as light as keV 
or a few GeV, as some experiments suggest, this scenario is ruled out. However for DM mass 
around 30 GeV, allowing small deviations from direct detection, the HTM with vectorlike 
leptons provides a viable explanation. This analysis assumed the DM candidate to be light 
and set an upper bound of 108 GeV, by the choice of the mass of the lightest verctorlikc 
charged lepton. One can extend this scenario to a more complicated one, involving several 
DM particles. This model would less constrained, but it looses the predictability of the 
simple scenario presented here. Given the importance of DM in understanding the universe, 
and the effort going into direct and indirect detection, and into collider experiments, simple 
models such as this one can help elucidating the nature of DM. 
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